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1. Introduction 

Various surface defects have been observed and investigated by scanning tunneling 
microscopy (STM) [U |2l [3l H]. The interference of the surface electron waves results 
in an oscillatory dependence of the tunneling conductance measured as a function of the 
separation between the STM tip and the defect. Remarkable manifestations of quantum 
interference were observed in artificial structures built from single atoms on a clean 
metal surface, the so-called quantum corrals [5J. Magnetic adatoms on non- magnetic 
host metal surfaces are of special interest as they produce a characteristic many-body 
resonance structure in the differential conductance near zero voltage bias attributed to 
the Kondo effect [H EJ [TJ [9] . The surface electrons waves contain the information of the 
magnetic impurity and by focussing the waves it has been possible to create a mirage 
image of the impurity [ID] ■ The shape of the resonance in the differential conductance, 
dl/dV, is usually asymmetric and is described by a Fano line shape [HI [121 US, HI]. 

In principle STM spectroscopy should also provide access to information on the 
structure of the metal below the surface. This possibility is based on the influence on 
the conductance caused by quantum interference of electron waves that are scattered 
by defects and reflected back by the contact. This effect was explored by Schmidt 
and coworkers ^5] for investigating subsurface bubbles of implanted gas in Al. The 
observation of interference patterns due to electron scattering by Co impurities in the 
interior of a Cu sample was reported by Quaas et al. [16]. Theoretically, the influence 
of single defects in the bulk of a metal on the quantum conductance of tunnel point- 
contact has been discussed in Refs. [TTl [18], [19]. In these papers it has been shown 
that the location of defects below the surface can be identified from the interference 
pattern in const ant- current STM images combined with the information obtained from 
the dependence of the conductance on the applied voltage. In the previous work of 
Refs. [TTl [iHl [IE] the scattering of electrons with a defect has been taken into account 
in the framework of perturbation theory. Such an approximation is valid as long as 
the strength of the electron - impurity scattering interaction is small. In the case of a 
magnetic defect at low temperatures {T <^Tk, where Tk is the Kondo temperature) the 
Kondo resonance results in a dramatic enhancement of the effective electron-impurity 
interaction [20] and the perturbation method becomes inapplicable. 

In this paper we present the quantum conductance G of the tunnel point contact 
in the vicinity of which a single point-like defect is situated, for arbitrary values of the 
scattering potential. We express the conductance by the means of a s-wave scattering 
phase shift 5q. The results describe the influence to the conductance of multiple 
scattering of the electrons by a single defect. Multiple scattering needs to be included 
even for a single defect because of electron reflection by the metal surface. This results 
in the appearance of harmonics in the dependence of G on the applied voltage and 
on the distance between the contact and the defect. We apply the analysis of the 
non-monotonic voltage dependence of the conductance specifically for the interesting 
problem of Kondo scattering, using an appropriate phase shift [21]. To our knowledge. 
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Figure 1. (a, 6) Model of the contact and (c, rf) illustration of the occupied energy 
bands in the two metal half-spaces for both signs of the applied bias eV . In panels 
(a, 6) the defect is placed at the point rg = (poj^o)- Electron trajectories are shown 
schematically. Note that we take the z-axis pointing downward. 

observation of subsurface Kondo impurities have not yet been reported in experiments, 
and the present analysis may guide future experimental investigations. 

2. Model and basic equations 

In our model of the system we represent the contact by an orifice of radius a centered 
at the location of the 'STM tip', r = 0. The orifice provides a tunneling window in 
otherwise impenetrable infinitely thin interface at 2; = between two metal half-spaces 
(Fig. [1]). The potential barrier at the plane of interface, z = 0, is taken to be described 
by a delta function, U (r) = Uof (p) 5 {z) , where p is the length of the radius vector p 
in the plane z = 0. The function / (p) ^00 in all points of the plane except in the 
contact, where / (p) = 1. At the point yq a defect described by the potential D (|r — ro|) 
is placed. 

We consider an almost ballistic configuration (the electrons are elastic scattered 
by the single defect only) and neglect electron-phonon scattering assuming the electron 
mean free path to be much large than the distance between the contact and the defect. 
In Ref. [22] the authors reported the observation of conductance oscillations at a voltage 
range up to l.beV at a temperature of 4.2K. Large bias voltages can be applied to small 
tunnel junctions created by STM or break-junction methods without significant heating 
of the electrodes. Because of the high resistance of the contact the current density 
remains small. Below we restrict our plots by the range eV < Sp. 

For the host metal we will consider a free electron model with an electron effective 
mass m* and a dispersion relation = ^^A;^/2m*, where k, and ey_ are the electron 
wave vector and electron energy, respectively. The electron wave function ■j/'k satisfies 
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the Schrodinger equation 



V^V^k (r) + [£k-t/(r)- V(r)]^k (r) 



D(|r-ro|)^k(r), 

where the V (r) is the apphed electrostatic potential. The function T/^k (p, z) satisfies 
boundaries conditions of continuity and of the jump of its derivative at the boundary 
2; = 0. We will assume that the transmission amplitude of electrons through the barrier 
in the orifice is small, 

t{k)^h^k/im*UQ]\t\<^l. (2) 

For small transparency t the voltage drop due to the applied bias is entirely localized 
at the barrier. The electric potential can be described by a step function, V {z) = 
VQ {—z) . As a result, the occupied energy bands in the half-spaces z > and z < 
are shifted by eV. We take the zero of energy, e = 0, to coincide with the bottom of 
the lower of the two bands, i.e. e = lies at the bottom of the band in the half-space 
2; > when eV > and at the bottom of the band in the half-space z < for eV < 0. 
At zero temperature electrons tunnel to the lower half-space (Fig. [Hc,d)) when eV > 0, 
and for eV < electrons can tunnel only to available states in the upper half-space. 

As shown in Refs. |23l [IT] Eq. ([T]) can be solved for arbitrary form of the function 
/ (p) in the limit |t| 0. To first approximation in the small parameter |t| -C 1 ([2]) the 
wave function tp^ (r) can be written as: 

(r) = ^ko (r) + tpki (r) , (3) 

where ip^i (r) ~ l/f^o- This latter part of the wave function ([3]) describes the electron 
tunnelling through the barrier and determines the electrical current. The first term 
in the Eq. ([3]) is the solution of the Schrodinger equation for the metallic half-spaces 
without the contact. It satisfies the boundary condition ■j/'ko (P; 0) = at the interface. 

For |t| <^ 1 the boundary condition for the jump of the derivative of the total wave 
function is reduced to 



-p ^kO 



z==FO 



^'"Vo/(p)4t)(p,0), (4) 



dz 

where V'kt^ ~ 0' -'-) wave functions for z ^ 0, k is the electron wave vector 

for electrons arriving in one half-space from the another half-space through the orifice, 
(|k| = ^k^ -2m*\eV\/h^). 

Thus, the function ipi^i (r) can be expressed by means of the solution ■i/'ko (r) . By 
using the Fourier transform of the wave function ([3]) we find 



(5) 



dp - 



{2nf J '^fip) 
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where k'^ = V/c^ — The electron wave function, ■z/'k (r) , which takes into account 
the scattering by the defect, can be expressed by means of the retarded Green function 
Gq (r', r; e) of the homogeneous equation at D = and [/ — oo. To first 
approximation in the transmission amphtude t ([2]) the integral equation for ■^/'ks (r) 
is given by 

V'ks (r) = (r) + 

2m* /■ , , , _ , „ (6) 



J dr'D{\r'-ro\)G+ (r, r'; e) V^ks (r 



where 



Gt (r, r'; e) = {/^« (A: |r - r'|) - (fc |r - , (7) 

r' = (p', —2;') . In Eq. ([7]) and below h^^^ (x) are the spherical Hankel functions. The first 
term in the braces is the Green function for free electrons in the infinite space and the 
second one takes into account the specular electron reflection from the interface. The 
functions tp^] (r) are the wave functions to zeroth and first order in t in the absence of 
the defect {D = 0). The electron wave function in the metal half-spaces is 

(r) = e**^" (e^'^^l"! - e"*'^^!"!) , (8) 

where k and kz are the components of the vector k parallel and perpendicular to the 
interface. The wave function ijj^^ (r) of the electrons that are transmitted through the 
contact has been obtained in Ref. [23]. In order to simplify further calculations we 



consider a point contact, taking the limit a 0. The solution ijj^^ (r) in this limit is 
given in Ref. [H] for any arbitrary anisotropic quadratic electron dispersion law ^k- For 
an isotropic band ^k = h^k'^ /2m* this takes the form, 

4?W = *(^-)^^^^"S"(*r)^ (9) 

Here, (r, 6', yj) are the spherical coordinates of the vector r, with 9 the angle between 
the vector r and the contact axis, k^ is the z— component of the vector k . The plane 
wave ([8]) is transformed into a spherical p-wave h'^p (kr) after scattering by the point 
contact. 

This model allows us to solve the three dimensional Schrodinger equation in the 
limit of small transparency of the barrier and find the analytical formulas for the 
conductance. Our method is similar to the widely employed tunneling Hamiltonian 
approach, where in the limit of small transparency of the barrier the distribution 
functions of electrons in the electrodes can be taken to be in equilibrium (Fermi 
functions) with chemical potentials shifted by the bias eV. For a barrier of finite width 
the electric field distribution changes which infiuences the nonlinear dependence of the 
conductance. This dependence becomes very important if the bias is comparable with 
the work function of the metal. For any three dimensional models of the potential 
barrier the dependence G {V) may be calculated only numerically. In our paper we 
did not make it our aim to investigate the intrinsic conductance of the tunnel junction 
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Go (y). The purpose of the work is to investigate the oscillatory and resonance additions 
to the conductance Gq {V) in the presence of a defect in the bulk of the metal, where 
the distribution functions are in equilibrium (in leading approximation in the barrier 
transparency). We believe that the part of the conductance related to the defect, which 
will be obtained in next sections, is correct, if the bias eV is less than Fermi energy sp- 



3. Scattered wave function in s-wave approximation 

Let D (|r — rol) be a spherically symmetric scattering potential which is finite in the 
point r = ro and tends to zero at a distance r/j <^ ro that is of the order of the Fermi 
wave length Ap- As is well known, s-wave scattering is dominant for scattering by a short 
range potential [21]. In order to express the wave function ([6]) by the s-wave phase shift 
5q we use the 'sharpness' of the function D (|r' — ro|), which essentially differs from zero 
only in a small region of the radius rjj near the point r' = ro. The main contribution to 
the integral in Eq. ([6]) comes from this region and the 'smooth' functions ip-^^s (r') and 

(A;|r-F|) can be taken outside the integral at the point r' = yq. For |r — ro| ^ 
the solution of Eq. IQ takes the form jl8], 



where 



^k. (r) ^ C (r) + (k) (ro) G+ (r, r^; e) , (10) 



Tik) = — ^ (11) 

Y{k)-ghi;^ {2kzo) 



_j_ m*ik 



Y{k)= J dr'D (r') (^^0 ,9 = j dr'D (r') . (12) 

Let us compare the wave function flTUl) with the formal solution ip^ (r) of the 
scattering problem for the spherically symmetrical potential D (|r — ro|) in infinite space 

(r) ^ (r) - (k) (ro) {k ko - r|) , (13) 

where ip^ and ip'^ are incident and scattered waves, and 

To (k) = 3^—-, (14) 

is the T matrix. Taking into account the relation between Tq and the s-wave phase 
shift 6o (k) 

^* To = ^e^^osin^o, (15) 



2TTh^ ^ k 
we rewrite the Eq. f fTTj) in the form 

T {k) = jr. . 16 

m*^k 1 + 1 (e^iSo _ 1) hW {2kz^) 

Note that the effective T-matrix flTBl) is an oscillatory function of the distance zq between 
the defect and the interface that results from repeated electron scattering by the defect 
after its reflections from the interface. 
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For a calculation of the current we should know the wave functions of the electrons 
transmitted through the contact, from one half-space to the other. For 2; > and 
eV > (i.e. for electron tunneling into the half-space in which the defect is situated) 
we find 

V'it^ (r) = <^ (r) - 

ik) (ro) {h^i' {k |r - rol) - {k |r -?o|)} • ^'^^ 

For 2; < and eV < (i.e. for electron tunneling from the half-space in which the defect 
is situated) the ip^^ (r) is written as 

(r) = <^ (r) + (ro) /^f ^ (kr) h['\kro). (18) 

Here V'ko 1 (r) ^"^^ T {k) are given by Eqs. ([H]), (E]) and (|T6l) . The wave functions ( fTTI) and 
ffTSjl have a completely different form: In the lower half-space the wave function f|T7|) is 
the superposition of the transmitted p-wave tp'^^ ~ /i^^^ (fc'"), dl]), and two s- waves, one of 
which, /ig^"* (A; |r — rol), is the wave scattered by the defect and other one hl^^ (fc |r — ro|) 
is the scattered wave, which undergoes reflection from the interface at 2; = (the wave 
moving from the 'image' defect placed in the mirror point Tq, |ro — FqI = 2zq). In 
the upper half-space there is only the p-wave ~ h^^^ (kr), the amplitude of which 
depends on the scattering on the defect because the wave incident to the contact is not 
a plane wave in this case. 

4. Total current and conductance 

The tunneling current I(V) = I^~^\V) — I^~\V) is the difference between two currents 
flowing through the contact in opposite directions. Each of them can be evaluated by 
means of the probability current density integrated over a half-sphere of arbitrary radius 
r, centered at the point contact r = and covering the contact from the appropriate 
side, and integrating over all directions of the electron wave vector. In this case the 
integrated probability current density jj^'^\v) is written as 

)(\/) = - — f dne (±z) [ dn^Q i±kz 

m* 




Im ^|.r(r)-^i^ , (19) 



where dfl and dfl^ are elements of solid angle in the real and momentum spaces, 
respectively. The total current through the contact is 

m = ^J dke \jt\v)f^ (.k - eV) X 

(2^) { ^ (20) 

(1 - /f (5k)) - Jt\v)h (ek) (1 - /f (^k - eV)) 
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where /p (^k) is the Fermi function. At zero temperature only one of the terms in square 
brackets in Eq. fl20l) differs from zero, i.e. only in one of the half-spaces states are 
available for tunneling, depending on the sign of the bias. Using the wave functions ( fTTI) 
and (fTSl) . after integration through Eq. (fT9l) the electrical current /*^^)(V^) at \eV\ < e-p 
and T = takes the form 



JW(\/) 



eha^ 
367rm* 



^yk^+2m■'\eV\/h? 



t{k) (1 + $(A;W)) 



(21) 



where the integration is carried out over the absolute value of the wave vector k 
within the interval \eV\ of allowed energies. We define fc^^-* = fc, /c*^"-* = k = 
^Jk^ — 2m* \eV\ /K^, kp is the Fermi wave vector. 



and 



$ (k) = sin5o-| [12ji (fcro) {-yi (kro) cos^o + 

{ji (^^^o) (jo C^kzo) -l)+yo i2kzo) yi (kro)} sin 6o) + 
6 (1 - Jo i^kzo)) (kro)-' (1 + (krof) sin^o] , 



D = 1 + 2 sin (5o X 



(22) 



jo{2kzo) sin^o - yoC^kzo) cos5o 



(23) 



,2(2A;zo)' 

and ji{x) and yi{x) are the spherical Bessel functions. From Eq. ( 12T|) it follows that 
the current-voltage dependence need not be symmetric in voltage in the presence of a 
defect. 

The differential conductance G = dl/dV for \eV\ < ep and for eV > is, given by 



G{V) = Go 
and for eV < 0, 

G{V) = Go 



kp 



q{V)(l + ^kp] 



kp 



dkk^^ {k) 



kp 



(24) 



q{V) + %{kp) 



kp 



kp 



dkk'k'^ (k) 



kp 



(25) 



Here kp = ^/kp + 2m*eV/h? and. 



qiV) 



Go=\t{kp)f 



h?kl 



2m* \eV\ 1 (2m* \eV 

(kpa)"^ 
36nh 



(26) 
(27) 
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is the conductance of the tunnel point-contact in the absence of a defect in the hmit 
V ^ At low voltage the conductance can be expressed as an expansion in the 
parameter 1/ {k-pZo) < 1, 

v2 



G(0) = Go i 1 + 12§^ (-1)" T^l^x 



n=l 



ill 



sin (2/cf (ro + («-!) 2:0) + nd^ 



cos (2A;f (ro + (n - 1) zq) + n^o) 



(28) 



The second term in the Eq. fl28|) corresponds to the sum over n scattering events by the 
defect and n — 1 reflections by the surface. If we keep only the term for n = 1 Eq. (l28ll 
is consistent with the results obtained by perturbation theory previously pTl [T8| [19]. 



5. Discussion and application to Kondo scattering 

The expansion fl28|) of the conductance G demonstrates that as a result of multiple 
scattering the conductance Go, Eq- ( pTl) . of the tunnel point contact becomes modified 
with oscillatory contributions AG„, which at 1/ (kpZo) <^ 1 and zq ~ ro is of order 

AG„ ~ ^, sin (2fcF (ro + {n - 1) zq) + n6o) , (29) 

(A^pro) 

where n = 1, 2... is the number of scattering events on the defect placed at a distance ro 
from the contact (and at a distance zo from the interface), and (n — 1) is the number of 
reflections by the interface. The argument of the sine function in Eq. fl29l) corresponds 
to the phase the electron accumulates while moving along a semi classical trajectory. 
In Fig. [I](a,b) such trajectories are illustrated for the case of scattering twice by the 
defect and one specular reflection by the interface. For eV > (Fig{T^) this trajectory 
consists of a segment (labelled 2) passing through the contact and arriving at the defect, 
two line segments (3 and 4) connecting the defect and the interface (these segments are 
perpendicular to the interface because only along such trajectory the electron can return 
to the defect and undergo the second scattering), and the part (5) from the defect to the 
contact. After specular reflection from the contact this wave interferes with the partial 
wave (1) that is directly transmitted through the contact. 

When eV < (Fig. [TJd) a wave incident to the contact (trajectory 2) is partially 
reflected from the contact. The electron moving along the trajectory 3 from the contact 
to the defect is partially scattered towards the interface (line segments 4) where it 
undergoes specular reflection from the interface (5) and comes back to the defect, from 
which it returns to the contact via trajectory 6. Tunnelling through the contact this 
partial wave interferes with the partial wave that is directly transmitted (1) in the half- 
space z < 0. At each scattering on the defect the electron acquires an additional phase 
shift Sq. The phase shift A0 between the two interfering partial waves for an electron 
with wave vector k is A0 = 2fcro + 2kzo + 26o- Because the maximum value of the 
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Figure 2. Voltage bias dependences of the normalized conductance corrections 
/S.G{V) / Go for a magnetic and a nonmagnetic impurity calculated from Eqs. (p4)) 
and ([25|l . We have used the parameters £k = O.Qep, ^bTk — O.Olep, = O.IAf/Ztt, 
po = 0, and zo = 20AF/27r. 



electron wave vector depends on the applied voltage the conductance oscillates as the 
function of eV . 

{The differential conductance, as the derivative of the current, discriminates a 
bound of the energy interval, which depends on the bias eV, i.e. for eV > the period 
of oscillations is defined by the energy ep + eV and for eV < - by the energy ep — \ eV\ . 

However, the current voltage characteristics is not symmetric relative to the point 
V = 0. This asymmetry results from the dependencies of the phase shift 6o{kp) and 
the absolute value of the wave vector k -p = \/kp + 2m*eV/ h? on the sign of eV . The 
physical origin of this asymmetry is that the scattering depends on the electron energy 
in the lower half-space, which is different for different directions of the current. 

The dependence doik) on k is defined by the form of the scattering potential U{r). 
To illustrate the obtained results for an s-wave phase shift we use the following model 
function 0121], 



The first term in Eq. (l30l) describes the resonant scattering on a Kondo impurity level 
£k (^k is the Kondo temperature). For Ek the effective electron scattering cross 

section acquires a maximum value corresponding to the Kondo phase shift 5ok = 7r/2 
[2U] . For a non-magnetic impurity this term is absent. The second term takes into 
account the usual potential scattering. For simplicity we use the s-wave phase shift for 
a hard sphere potential of radius {kprp, < 1). Inelastic scattering by the magnetic 
defect can be taken into account in the scattering formalism by introducing an imaginary 
part of the phase fl30|) . 

Figure [2] shows the dependences of the corrections to the normalized conductance 
AG(y)/Go = {G(y) — Go(y)) /Gq resulting from the scattering by a defect placed on 
the contact axis for a magnetic and a nonmagnetic impurity. The figure illustrates 




(30) 
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Figure 3. Difference 5Gk{V)/Gq between tlie voltage dependences of the conductance 
for a magnetic and a non-magnetic impurity. We have used the parameters Ek — O.Qep, 
/fcBTk = O.Olep, and td = 0.1AF/27r. 

the appearance of a Kondo anomaly in the conductance seen as an extremum in 
the differential conductance, GiV), near the bias cVk corresponding to the resonance 
condition + cVk — £:k = 0. The plots show a slowly increasing background on top of 
the oscillating AG'(y) dependence. The background arises from the integral terms in 
Eqs. (1241) . fl25l) . which take into account the contribution of all available states within 
interval \eV\. The monotonic part in AG(V^) is more pronounced in the case of Kondo 
scattering, which gives a large contribution to this part at any voltage. 

It is interesting to observe that the sign of the Kondo anomaly depends on the 
distance between the contact and the defect rg. This distance in combination with the 
value of the wave vector k-p determines the period of oscillation of AG'(V^), which is 
indeed a non-monotonic function of /cp'^o- If the bias eVx coincides with a maximum in 
the oscillatory part of conductance the sign of the Kondo anomaly is positive and vice 
versa, the negative sign of the Kondo anomaly is found at a minimum in the periodic 
variation of AG. 

In Figl3]we present the difference 5G}<^{y) jG^ = (AG^ — AG„)/Go between voltage 
dependences for a magnetic AGm and a non-magnetic AG„ impurity, having the same 
potential scattering strength. The plots in the Fig. [3] show the evolution of the shape 
of the Kondo anomaly for several values of the distance between the contact and the 
impurity, placed on the contact axis. The change of distance changes the periodicity 
of the normal-scattering oscillations which is illustrated to lead to a changing of sign 
in the Kondo signal. A similar dependence of the differential conductance with the 
distance between an STM tip and an adatom on the surface of a metal has been obtained 
theoretically in Refs. [131 125] in the terms of the Anderson impurity Hamiltonian |26j . 
Note that we obtained the Fano-like shape of the Kondo resonance in the framework a 
single-electron approximation while in Refs. flSl [25] the many-body effects were taken 
into account. 
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Figure 4. Comparison of the oscillatory parts of the conductance AG(l^)/Go 
calculated by using the Eqs. ([24l [25]) (full curves) and by means of results obtained in 
the framework of perturbation in the electron- impurity interaction (dashed curves), a 
- non-magnetic defect; h - magnetic defect. We have used the parameters ek = O.Oep, 
/cbTk = O.Olep, rj) = 0.1AF/27r, po = 0, and zq = 5AF/27r. 

Figure m illustrates the importance of multiple scattering for this problem. It shows 
the oscillatory parts of the conductance AG{V)/Go calculated by using the Eqs. (IMl 
l25l) in comparison to results obtained in the framework of perturbation in the electron- 
impurity interaction [171 HE], i.e. neglecting multiple electron scattering. While for the 
non- magnetic impurity ((Fig. H^)) the difference between two curves is small it is seen 
that for a magnetic impurity (Fig. |Dd) the perturbation method does not describe the 
conductance correctly in a region of the Kondo resonance. For nonmagnetic impurities 
multiple scattering has a negligible effect due to the smallness of contributions of the 
multiple scattering paths described by the the parameter {k f^o)~^ , which is no longer 
true near the Kondo resonance, where the increasing of the scattering amplitude is the 
dominant effect. 

6. Conclusion 

We have studied the influence of multiple electron scattering by a single defect on the 
current through a tunnel point- contact. In the approximation of s-wave scattering by 
the defect a general expression for the conductance G has been found ( !24l) . (!25ll . The 
results obtained have been analyzed for the model s-wave phase shift ( pOj) describing 
the Kondo scattering by a magnetic impurity. We demonstrated that taking multiple 
scattering into account is most essential near voltage values corresponding to the Kondo 
resonance condition + eV = Ek- It is found that the the shape as well as the sign of 
the Kondo anomaly depends on the position of the defect. This dependence results from 
quantum interference of partial waves directly transmitted through the contact with the 
partial wave scattered by the defect and reflected by the interface. The phase shift 
between the two waves produces the oscillations of the conductance. A maximum in the 
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regular oscillation of G leads to a positive sign of the Kondo anomaly at that position, 
while a minimum produces a negative sign. These results may be exploited in future 
experiments for detecting and investigating the Kondo effect of individual impurities in 
the bulk of a host metal. 

Acknowledgment s 

Ye. S. A. is supported by the INTAS grant for Young Scientists (No 04-83-3750) and Yu. 

A. K. was supported by a NWO visitor's grant. This research was supported supported 

partly by the program "Nanosystems nanomaterials, and nanotechnology" of national 

Academy of Sciences of Ukraine. 
[1] Crommie M.F., Lutz CP. and Eigler D.M., 1993, Science, 262, 218. 

[2] Sprunger P.T., Petersen L., Plummer E.W., Lagsgaard E. and Besenbacher F., 1997, Science, 275, 
1764. 

[3] Kolesnychenko O.Yu., de Kort R., Katsnelson M.I., Lichtenstein A.I. and van Kempen H., 2002, 
Nature, 415, 507. 

[4] Hofmann Ph. , Brincr E.G., Doering M., Rust H.-P., Plummer E.W. and Bradshaw A.M., 2002, 

Phys. Rev. Lett, 65, 115420. 
[5] Fiete G.A., Heller E.J., 2003, Rev. Mod. Phys., 75, 933. 

[6] Madhavan V., Chen W., Jamneala T., Crommie M.F. and Wingreen N.S., 1998, Science, 280, 
567. 

[7] Li J., Schneider W.-D., Berndt R. and Delley B., 1998, Phys. Rev. Lett., 80, 2893. 

[8] Chen W., Jamneala T., Madhavan V. and Crommie M. F., 1999, Phys. Rev. B, 60, R8529. 

[9] Knorr N.,Schneider M.A., Dickhoncr L., Wahl P. and Kern K., 2002, Phys. Rev. Lett., 88, 096804. 
[10] Manoharan H.P., Lutz CP. and Eigler D.M., 2000, Nature, 403, 512. 
[11] Fane U., 1961, Phys. Rev., 124, 1866. 

[12] Plihal M. and Gadzuk J.W., 2001, Phys. Rev. B, 63, 085404. 

[13] IJjsaghy O., Kroha J., Szunyogh L. and Zawadowski A., 2000, Phys. Rev. Lett., 85, 2557. 
[14] Wahl P. et.al., 2004, Phys. Rev. Lett., 93, 176603. 

[15] Schmid M., Hebenstreit W., Varga P. and Crampin S., 1996, Phys. Rev. Lett, 76, 2298. 

[16] Quaas N., Wenderoth M., Weismann A., Ulbrich R.G. and Schonhammer K.,2004, Phys. Rev. B, 

69, 201103(R). 

[17] Avotina Ye.S., Kolesnichenko Yu.A., Omelyanchouk A.N., Otte A.F. and van Ruitenbeek J.M., 

2005, Phys. Rev. B, 71, 115430. 
[18] Avotina Ye.S., Kolesnichenko Yu.A., Otte A.F. and van Ruitenbeek J.M., 2006, Phys. Rev. B, 74, 

085411. 

[19] Avotina Ye.S., Kolesnichenko Yu.A., Otte A.F. and van Ruitenbeek J.M., 2007, Phys. Rev. B, 75, 
125411. 

[20] Abrikosov A. A., Fundamentals of the theory o J metals, North Holland, 1988. 
[21] Schneider M.A., Vitah L., Knorr N. and Kern K., 2002, Phys. Rev. B, 65, 121406. 
[22] C. Untiedt, G. Rubio Bollinger, S. Vieira, and N. Agrai't, 2000, Phys. Rev. B, 62, 9962. 
[23] Kulik I.O., Mitsai Yu.N. and Omelyanchouk A.N., 1974, Zh. Exp. Tear. Fiz., 63, 1051. 

[24] Davidov A.S., Quantum Mechanics, Pcrgamon Press, 1976. 

[25] Chiung-Yuan Lin, Castro Neto A.H. and Jones B.A., 2005, Phys. Rev. B, 71, 035417. 
[26] Anderson P.W., 1961, Phys. Rev., 124, 41. 



